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Abstract — Partial decoding has the potential to achieve a larger 
capacity region than full decoding in two-way relay (TWR) 
channels. Existing partial decoding realizations are however 
designed for Gaussian channels and with a static physical layer 
network coding (PLNC). In this paper, we propose a new solution 
for joint network coding and channel decoding at the relay, 
called pairwise check decoding (PCD), for low-density parity- 
check (LDPC) coded TWR system over block fading channels. 
The main idea is to form a check relationship table (check- 
relation-tab) for the superimposed LDPC coded packet pair in 
the multiple access (MA) phase in conjunction with an adaptive 
PLNC mapping in the broadcast (BC) phase. Using PCD, we then 
present a partial decoding method, two-stage closest-neighbor 
clustering with PCD (TS-CNC-PCD), with the aim of minimizing 
the worst pairwise error probability. Moreover, we propose the 
minimum correlation optimization (MCO) for selecting the better 
check-relation-tabs. Simulation results confirm that the proposed 
TS-CNC-PCD offers a sizable gain over the conventional XOR 
with belief propagation (BP) in fading channels. 

Index Terms — Two-way relaying, block fading channel, LDPC, 
physical layer network coding, partial decoding, pairwise check 
decoding, closest-neighbor clustering. 



L Introduction 

Two-way relaying, where two source nodes exchange infor- 
mation with the help of a relay node, has recently gained a 
lot of research interests |3|-[12J. It is shown able to overcome 
the half-duplex constraint and significantly improve the system 
spectral efficiency in relay-based cooperative networks. Upon 
receiving the bidirectional information flows, the relay node 
combines them together and then broadcasts to the two desired 
destinations. The operation at the relay resembles network 
coding [13 K a technique originally developed for wireline 
networks. It is thus often referred to as physical layer network 
coding (PLNC) HJ or analog network coding (ANC) 18|. 

Among the various two-way relaying strategies, the two 
practical and efficient ones are known as amplify-and-forward 
(AF) and decode-and-forward (DF), similar to those in one- 
way relaying. Different from one-way relaying, the DF strat- 
egy for two-way relaying further includes full DF O, ifTSl 
and partial DF II16II - II18L This is because the combining 
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process operated at the relay is a many-to-one mapping (e.g. 
I0l = 000 = 0). As a result, it is not necessary for the 
relay to fully decode the message pair before combining them 
together. Being a simple realization of partial DF, the denoise- 
and-forward (DNF) strategy proposed in |19| demonstrates 
significant performance gain over the full DF. It is also known 
that partial decoding has the potential to achieve a much larger 
rate region than full decoding, even though its capacity region 
is still unknown. Consequently, it remains as a fundamental 
and challenging task to realize the potential of partial decoding 
in two-way relay (TWR) channels through practical coding 
and modulation techniques. 

Several works have been reported on the implementations of 
partial DF for channel-coded TWR systems, and they are also 
known as joint network-coding and channel-coding (JNCD) 
design in the literature. An intuitive method is to utilize the 
fact that the network-coded (eg. XOR or modulo addition) 
codeword pair is also a valid codeword given the same linear 
code (e.g. LDPC or Lattice codes) applied at both source 
codes ||20 | -||23 l . In this method, the relay first computes the 
probability of based on the received superimposed 
signal yc during the multiple-access (MA) phase, where xa 
and xb are the symbols after channel coding and modulation 
from source A and source B, respectively, and then apply soft- 
decoding to decode the associated network-coded information 
symbol pair sa® sb- We refer to this method as partial DF 
based on conventional XOR. However, this method discards 
useful information related to the decoding of the whole packet 
sa®sb during the mapping from signal yc to the probability 
of Xa ® Xb- ^ more advanced partial decoding method is 
to exploit the Euclidean distance profile of the superimposed 
symbol pair after going through the noisy channel in the MA 
phase ll24l . The relay first decodes the arithmetic- sum of the 
coded symbol pair ca ^ cb and then map it to sa® sb for 
broadcasting. The authors in [24J show that this partial DF 
method based on arithmetic- sum provides higher decoding 
gain than the one based on conventional XOR. Note that 
both aforementioned methods are designed specifically for 
symmetric and Gaussian channels. 

For TWR channel with fading, the conventional XOR does 
not always work well due to the undesired phase and amplitude 
offset between the two TWR channels in the MA phase. 
Authors in 1 25 1 therefore proposed an adaptive network coding 
with respect to the instantaneous channel fading, named as 
closest-neighbor cluster (CNC) mapping. Compared to the 
conventional XOR, the CNC mapping obtains a higher end-to- 
end throughput. To further ensure reliable communication, the 
authors extended this method for convolutional-coded system 
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in l|26ll and discussed the code design based on trellis-coded 
modulation (TCM) f27, Section 8.2]. However, this TCM- 
based CNC mapping requires to change the coding structure at 
the two source nodes and switches two transmission protocols 
(CNC DNF and pseudo AF) in order to exploit the system 
performance. 

In this paper, we propose a new relay channel decoding 
solution, called pairwise check decoding (PCD), for LDPC 
coded TWR fading channels. The main idea is to form a check 
relationship table (check-relation-tab) for the coded symbol 
pair {ca-,cb) by taking both the employed LDPC codes and 
the adaptive PLNC mapping into accounts. The proposed 
PCD algorithm is universal for any adaptive PLNC mapping 
and does not require the LDPC codes at the two sources 
being identical. It offers a practical and efficient approach 
to realizing the promising DNF TWR strategy with advanced 
channel coding. With the aim of maximizing the minimum 
Euclidean distance (MED) between any two codewords, we 
present a partial decoding method, two- stage CNC with PCD 
(TS-CNC-PCD), for TWR fading channels. The proposed TS- 
CNC-PCD is appropriate for any choice of constellation size. 
Moreover, a kind of correlative rows optimization, named as 
the minimum correlation optimization (MCO), is proposed 
for selecting the better check-relation- tabs. Simulation results 
confirm that the proposed TS-CNC-PCD has significant coding 
gains over the conventional XOR with belief propagation (BP) 
decoding algorithm. 

The rest of the paper is organized as follows. In Section 
II, we present the channel coding model for TWR block 
fading channels. Section III analyzes the lower bound of the 
outage probability. The design criterion of adaptive codeword 
mapping is interpreted in Section IV. In Section V, we propose 
the PCD algorithm in detail. Section VI present a two-stage 
CNC mapping based on PCD decoding, named as TS-CNC- 
PCD, in terms of optimizing the MED. The convergence 
behaviors and the coding gains of the proposed TS-CNC-PCD 
method are simulated in Section VII. Finally, we conclude the 
paper in Section VIII. 

II. Channel coding model for two-way relay 

FADING CHANNELS 

We consider a TWR fading channel where two source 
nodes, denoted as A and B, exchange information with the 
help of a relay node, denoted as C. We assume that all 
the nodes operate in the half-duplex mode. The channel on 
each communication link is assumed to be corrupted with 
block fading and additive white Gaussian noise (AWGN). For 
simplicity, we also assume the channel gains are reciprocal 
and unchanged during a whole packet transmission. 

The proposed channel coding model is shown in Fig. [T] 
where the communication takes place in two phases. First, 
the information packet from each source, denoted as S^, for 
i G {A^B}, is encoded individually by a traditional LDPC 
code with parity check matrix F^. Unlike the existing work, 
we do not impose the constraint that Ta and Tb must be 
identical. Instead, we only require that they have the same 
size and the same location of non-zero elements. We further 



assume that the encoder is operated in GF(g), where q G 
{2^, 2^, 2^,...}. Note that g-ary (g > 2) coding can improve 
the performance compared with binary coding. More details 
can be found in ll28l - ll3Qll and references therein. The encoded 
packet, Ci, is modulated by using g-ary modulation, denoted 
as Qq, such as g'-PSK or g-QAM, generating X^, and then 
transmitted simultaneously to the relay node. The n-th symbol 
of each packet is denoted as C>i{n) G Zq, = {0, 1, . . . , g — 
1}, and Xi(n) G Qg, respectively. The superimposed packet 
received by the relay, denoted as Yc is given by 

Yc = Hac^a + Hbc^b + Wc, (1) 

where Hu^ denotes the complex- valued channel coefficient of 
link from node i to node i\ and W^/ denotes complex AWGN 
with variance af, of node ^^ Therein, z, G {A, B, C}. 

We assume perfect symbol synchronization at the two 
sources and perfect channel estimation at the relay. After 
receiving the superimposed packet, the relay first computes 
the soft information (e.g. likelihood value) about the codeword 
after PLNC mapping, denoted as Cc = M.{Ca^Cb), then 
applies the proposed PCD to obtain the hard-decision (the 
details will be presented in Sections V and VI). Here, M. 
denotes a kind of adaptive PLNC mapping, which is designed 
on the codeword pair, and Cc(^) G Zq', where is the 
cardinality of the PLNC mapped symbol Cc(^), which should 
he q < q' < q'^. Note that no full decoding of Ca and is 
needed as an intermediate step. No extra channel encoding 
at the relay is needed either. Then, the relay broadcasts 
the modulated coded symbols of Cc, denoted as Xc, and 
mapping rule M{-) to two source nodes. The received signals 
at the nodes A and B are respectively written as 

Ya = Hca^c^Wa; ... 
Yb = Hcb^c^Wb. ^ ^ 

Each source node, say B{A), computes the likelihood of 
the desired information Ca{Cb) from the received symbols 
Yb{Ya) by using the inverse PLNC mapping rule with the 
help of its self-information Cb{Ca)- Lastly, the traditional 
LDPC decoding algorithm, e.g. BP, is applied, the output of 
which is the desired information packet Sa(Sb). Note that 
each source node should know the check matrix of the other 
source. 

III. Analysis of outage probability 

In this section, we derive the system outage probability of 
TWR fading channel, which serves as a good approximation of 
the achievable frame error rate (PER) in the limit of infinite 
block length [31 1. Here, the system is said to be in outage 
if the achievable sum-rate falls below a target. Since the 
capacity region of two-way relaying with partial decoding is 
still unknown |[T6ll - |[T8ll . we resort to the capacity outer bound 
as follows mi Theorem 2], based on which a lower bound of 
the outage probability can be obtained. 

A : \Rab < min {^Cac. (1 - P)Ccb). 

^ ^ (3) 

Rba < min {/3Cbc. (1 - P)Cca)\. 
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Source (A) / Destination (B) 



Relay(C) 



Source (B) / Destination (A) 




Fig. 1: Channel coding model for TWR block fading channels. 



where P is the time sharing parameter, Rij and Cij denote the 
instantaneous data rate and channel capacity of the link from 
node i to node j, for z, j G {A^B^C}, respectively. 

After simple manipulation of the constraints in (|3]Q, we 
obtain the following linear inequalities about Rab and Rba 



Rab 
Cac 
Rab 
Cac 
Rab 

QCB 

Rbc 
Cbc 



Rab 

CcB 

Rba 

CcA 

Rba 
Cbc 

^BC 
CcA 



<1, 
<1, 

< 1. 



(4) 



Let us further assume that the TWR channels considered 
here are reciprocal, i.e. Cij = Cji for i^j e {A, 5, C}. From 
(O, we can easily obtain the upper bound of the maximum 
sum-rate for the considered TWR channels 



max Rab + Rba = min(CAc, Cbc), (5) 

iRAB,RBA)eA 



which is also given in |19|. Therein, each of the terms Cij, 
i^j e {A, 5, C}, is the channel capacity of a traditional point- 
to-point channel with input alphabet Xij G Qq and received 
signal yij = aijXij + Wij, where Wij ~ A/'(0,<j^) and aij 
denotes a real- or complex- valued channel coefficient of the 
link from node i to j with Ejla^jp} = 1. 

With the further assumption of equiprobable channel in- 
puts, extending the well-known formula for the capacity of 
continuous-valued Gaussian channels |33, Eqs. 3-5] to the case 
of block fading channels yields 



Q,K-) = iog2(^)-i E E 



log2 E exp 

n=0 



m=0 

\yij-aijx'^j\'^-\yij-aijx'^\'^ 
2^2 



(6) 



in bit/channel use. Here, E represents expectation over yij 
given Xij = x^j and aij, with x'^ being an element of the 

modulated signal sets | Qq : x^j , xjj , . . . , x^^"^ | . 

In addition, we denote the target sum-rate of overall system 
as Sr- In our considered channel coded TWR model, given 

^The detailed derivation is similar to that described in (32). 



that the same LDPC code rate, denoted as r, and the same 
constellation are employed by the two source nodes, we have 

Sr = rlog^{q). (7) 

Then the outage probability can be lower bounded as 

Pout 

= p(^mmi^CAc{(^Ac),CBc{oiBc)^ <^r), 

which can be easily evaluated by Monte Carlo averaging over 
the block fading coefficients and the AWGN. 

IV. Design CRITERION OF A^(Ca, Cb) 
Similar to the uncoded system in f25\, a necessary condition 
for successful decoding at two sources/destinations in coded 
TWR system is for the adaptive PLNC mapping to satisfy the 
exclusive law: 

M{Ca, Cb) + M{C'a. Cb), for ^\\{Ca ^ C'j,, C^}, 
M{Ca. Cb) + M{Ca. Gj,)M all{CA, ^ e^}- 

(9) 

Given the above necessary condition, we next discuss the 
design criterion of A1(Ca, Cb) to optimize the system error 
performance. We aim at minimizing the pairwise error prob- 
ability (PEP) between two distinct codewords and Cj^ 
in the MA phase because the MA interference dominates the 
whole system performance. Therein, Cj^ and Cj^ denote the 
corresponding codewords generated by the l^^ and V^^ code- 
word pairs (Ca,Cb) with certain selected PLNC mapping 
A^(-) respectively. Note that the PEP considered in this work 
is defined over codewords, whereas ll25l was concerned with 
the PEP over uncoded symbols. For TWR fading channels, the 
PEP between C^ and Cj^ conditioned on the instantaneous 
channel gain pair {Hac^ Hbc} is given by Oi, p. 265] ll27l 
Section 5.5] |35, Eq. 7] EH Eqs. 5-6] 

Pe{C'c^^C I {HaC.Hbc}) 

= Pe [m{C^a^ C^) ^ M(C^, eg I {Hac. Hbc}) 

'C,cDI,\{Hac.Hbc}\ 
2^^ )' 

(10) 

where and are the energy per coded symbol and the 
variance of Gaussian noise at relay node C, respectively. 
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Q{-) is the Q-function. Df^, represents the squared EucHdean 
distance between the codewords Cj^ and Cj^ and is the 
function of two channel coefficients {Hac-,Hbc} given by 



D. 



N 

E 

n=l 



(11) 



\xi{n) - xi> {n) . 
In (fTTI) , N denotes the length of transmitted codewords and 

(12) 



xi(n) = HAcx\(ri) + Hbcx^b^^)^ 
xi'{n) = HAcXAin) + Hbcx^b^^)- 



From ([TOl) , it is clear that to minimize the pairwise er- 
ror probability, one should design a mapping rule Cc = 
A^(Ca,Cb) that can maximize the MED between any pair 

C' ^c) 



of codewords (Cj^, C[^). We denote it as 



mm 



(13) 



Note that the optimal strategy of the mapping rule is adaptive 
with respect to the channel conditions. 

V. Pairwise check decoding (PCD) 

After introducing the design criterion of the adaptive PLNC 
mapping M{-) for coded TWR systems in the previous sec- 
tion, we present a general relay decoding framework, named as 
pairwise check decoding, for any given M{-) in this section. 

A. Check relationship table (check-relation-tab) at relay 

It is clear that if the conventional XOR mapping is applied, 
we can easily construct a virtual parity check matrix Tc for 
the codeword Cc based on Ta and at the relay and then 
decode Cc using the traditional BP algorithm |36|, |37|. It 
is because the referred XOR mapping is linear. However, if 
an adaptive PLNC mapping is applied, which is much more 
likely to be non-linear, such virtual parity check matrix Tc 
in explicit form cannot be found. Instead, one has to resort 
to the constraint relationship regarding the codeword pair. In 
this subsection, we introduce a so-called check-relation-tab to 
describe such codeword pair constraints. 

At first, we set forth some notations. We assume that two 
LDPC codes Ta and on GT{q) are used at the two 
sources. The two parity check matrices have the same size 
of M X and the same locations of non-zero's, where N 
is the codeword length and M is the number of parity check 
symbols. Let and F^^ denote the elements at the m- 
th row and n-th column of F^ and F b ■> respectively. = 
{n : Frnn 7^ 0} dcnotcs the set of column locations of the 
non-zero's in the m-th row; M^\^ = {n' : Vran' 7^ 
denotes the set of column locations of the non-zero's in the m- 
th row, excluding location n. Likewise, Nn = {m : Vmn 7^ 0} 
and A^n\m = {j^' • ^m'n 7^ dcnotcs the set of 

row locations of the non-zero's in the n-th column and those 
excluding location m. Note that in the following sections 
the arithmetic operations are all in GF(g') unless specified 
otherwise. 



1 ) Check-relation-tab: Considering the encoding character- 
istics of F^ and F^, we have the set of parity check equations, 
for each m = 1, . . . , M, as follows: 

0, 

(14) 



A 

mn 



E c^(n)xr 

E CB(n)xr^„ 



0. 



For each set of the above check equations, we construct 
one check-relation-tab. The m-th check-relation-tab essen- 
tially characterizes the joint constraint of the symbol pairs 
{C^(n), 0^(71)} at locations n G M^. It consists of two 
parts, one for virtual encoder, and the other for PCD decoder. 

In virtual encoder, without loss of generality, we assume 
that the symbol pair {Ca(^), Cb(^)} at a random location 
n is an unknown parity check symbol pair while those 
{C^(n'), Cb(^0} ^1 other non-zero locations n' G M^\^ 
are the known information symbol pairs. The parity check 
symbol pair and the information symbol pair are assumed 
as the symbol pairs composed of the unknown parity check 
symbols and the known information symbols respectively. 
We obtain the possible values for {CA(n), 0^(77,)} at the 
given n base on (O through enumerating all values for 
{C>a{'^')-,^b{'^')} at locations n' G M^\^. Then, the sym- 
bol pairs {Ca(^), Cb(^)} at all locations n G Mm are 
mapped to Cc(^), according to a given PLNC mapping 
rule Cc = A^(Ca,Cb). Since the number of symbol pairs 
{CA(n), Cb(^)} mapped to each element Cc(^) may not 
be the same, the probability of occurrence for each element 
Cc(^) should be computed separately. 

For PCD decoder, without loss of generality, we assume 
that the element Cc(^) at location n is known. That is to 
say, the set of symbol pairs mapped to Cc(^) is given. 
Then, we should compute the probability of occurrence for 
the corresponding possible values Cc(^0 ^1 locations n' G 
Mm\n^ named as weighted factor Fw, by classifying the 
aforementioned probability of each element which is generated 
by the virtual encoder. 

In the following we shall use a toy example to demon- 
strate the detailed construction of the check-relation-tab as 
mentioned above. Consider the two source LDPC codes on 
GF(2) with 



Ta 



10 10 1 
110 10 
1110 
10 11 



where the code length is 6 and the row weight and column 
weight are 3 and 2 respectively. 

In order to indicate that our proposed PCD approach can 
deal with any PLNC mapping under the exclusive law ©, a 
special non-linear clustering of the codeword pair, denoted as 
Mnu is considered here and characterized as: 

{c^(n),CB(n)}^{Cc(n)}: 

{(0, 1)} ^ {a}, {(0, 0), (1, 1)} ^ W, {(1, 0)} ^ {c}. 

Each symbol pair |CA(n), Cb(^)| is mapped to one 
of the three elements {a,6, c}. Here, to avoid confusion. 
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Segmental Tanner Graph 

Fig. 2: The Tanner graph of the check functions versus codeword 
pair at the relay. 



we let {a, 6, c} indicate the decoding symbols (they are not 
necessarily equal to broadcasted symbols) during the period 
of operating PCD algorithm. Such mapping Mni may be 
appropriate when the channel ratio satisfies Hbc/Hac = — 1 
and the BPSK modulation is applied at two source nodes. 

The corresponding Tanner graph of the virtual LDPC code 
observed at the relay is shown in Fig. |2l The solid circles and 
squares denote the check functions at each source node, while 
non- solid circles and squares denote the transmitted symbols 
of a code. In like manner, the solid and non-solid ellipses 
denote the check functions and the corresponding received 
symbol pairs. denotes the m-th check function of the 
LDPC code at the node s, where me [0, 3], s G {A, B, C}. 

Next, we will derive the check function of from 
and /(f using the segmental Tanner graph in Fig. [2l which 
is separated from the panoramic Tanner graph. 1^ denotes the 
symbol pair (C^(n), Cb(^)), for n G [0,5]. Then, the first 
set of constraint equations associated with the 0^^ row of Ta 
and can be written as: 

C^(0) + C^(2) + C^(5) = 0, 

Cb(0) + Cb(2) + Cb(5) =0. ^'^^ 

Given the Mni mapping, the so-called check-relation-tabs 
for virtual encoder and PCD decoder are generated, as shown 
in Table HI For virtual encoder, without loss of generality, we 
assume that the symbol pairs Iq and I2 have been known, e.g. 
lo = a and I2 = c, as shown in Fig. [21 The possible value 
for the parity check symbol pair ^5 can be obtained by (fTSl) . 
It is important to note that each element constrained by the 
pairwise check functions corresponds to one symbol pair 
in {/^, n e {0, 2, 5}}, the range of which is in {a, 6, c}. 

Through reverse derivation, the check-relation- tab for PCD 
decoder are given. Therein, Fw is the weighted factor which 
should be multiplied with the probability of the two elements 
behind. Then, we just need to compute the probability of 



TABLE I: Check-relation-tabs of /q for binary LDPC codes with 
BPSK modulation 



Virtual Encoder 


PCD Decoder 


(0,0) 


(2,2) 


(5,5) 


T-p =rTr 

(5,5) 




(0,0) 


— 77^ — f^^^ — 
(2,2) 


a 


(0,1) 


a 


(0,1) 


lb 


(0,0) 


a 


0.5 


a 


b 


a 


(0,1) 


b 


(0,0) 


0.5a 


(0,1) 


a 


0.5 


b 


a 






(1,1) 


0.5c 


(1,0) 


a 


0.5 


b 


c 


a 


(0,1) 


c 


(1,0) 


lb 


(1,1) 


a 


0.5 


c 


b 


b 


(0,0) 


a 


(0,1) 


0.5a 


(0,1) 


b 


1 


a 


a 




(1,1) 




0.5c 


(1,0) 


b 


1 


a 


c 


b 


(0,0) 


b 


(0,0) 


lb 


(0,0) 


b 


1 


b 


b 




(1,1) 




(1,1) 


(1,1) 


b 


1 


c 


a 


b 


(0,0) 


c 


(1,0) 


0.5a 


(0,1) 


b 


1 


c 


c 




(1,1) 




0.5c 


(1,0) 


c 


0.5 


a 


b 


c 


(1,0) 


a 


(0,1) 


lb 


(1,1) 


c 


0.5 


b 


a 


c 


(1,0) 


b 


(0,0) 


0.5a 


(0,1) 


c 


0.5 


b 


c 






(1,1) 


0.5c 


(1,0) 


c 


0.5 


c 


b 


c 


(1,0) 


c 


(1,0) 


lb 


(0,0) 











occurrence of the corresponding possible values {a, 6, c} for 
lo and I2 one by one. Interestingly, using the binary LDPC 
codes or the g-ary {q > 2) codes with the non-zero elements 
{7^, ... ,7^ G Zq}, we have the same probability distribution 
of occurrence among {in^n G {0,2,5}} for a given PLNC 
mapping, like = f^,,m ^ m' . 



B. Pairwise check decoding (PCD) algorithm 

After obtaining the check-relation-tabs, the PCD algorithm 
can be readily carried out. Note that the locations of non- 
zeros in the g-ary LDPC codes used at two source nodes 
are still valid here. We only change the check functions of 
symbol pairs from {f^^f^} to by the derived check- 
relation-tab. Define = Pr{Cc{n) = k\Yc{n)jWmn)', 
^mn = ^r(/^satisfied|Cc(n) = k.tnm)- Let tnm denotes 
the messages to be passed from symbol node Cc{n) to check 
node ; Wmn denotes the messages to be passed from check 
node to symbol node Cc{n). Suppose that {m/^}, k G Zq^, 
denote the index- set of rows, which have the same target 
value in the check-relation-tab for PCD decoder. Therein, the 
element at (n,n) is generated as {k} and rritab means the 
index of row. 

For a given PLNC mapping M, we compute for each [m, n] 
that satisfies Tmn 7^ 0. 
L Initialization 

Compute the initial value of each and each tnm as: 

'k 

E Pr((C^(n),CB(n))|Yc(n)); 

(CA(n),CB(n)):Cc(n)=/c ^ ^ 




where Pr\^{CA{'n)jCB{'n))\Yc{n)j is the probability of 
(C^(n), Cb(7i)) given Yc(n) is received. Obviously, we 
have u^^^ = u^^^ even if mi ^ m2,mi,m2 G N^. For 



6 



example, given Mni mapping^, we have 



nm 
'h 



exp 
exp 
exp 
exp 



-(Ycl 


[n)- 


-{Hac-Hbc)?\ 








-(Ycl 


[n)- 


-{Hac+Hbc))^\ 






2al ) 


-(Ycl 


[n)- 


-{-Hac-Hbc)? 








-(Ycl 


[n)- 


-{-Hac+Hbc)? 


2al 



(17) 



2. First half round iteration: from symbol node Cc(^) to 
check node 



^mn E Fw{mtah) H ^n'm 



(18) 



where tn'm{'^tah) denotes that /c of {u^,^^k G is the 
designated elements at the m^a6-th row in the check-relation- 
tab for PCD decoder. For Mni mapping, we have 

'^05 — 0-^('^00'^20 + '^00'^20 + '^00'^20 + '^00'^2o)5 

'^05 — '^00'^20 + '^00'^20 + '^00'^20 + '^00'^20 + '^00'^205 

'^05 ~ 0-5('^00'^20 + '^00'^20 + '^00'^20 + '^00'^2o)- 

(19) 

3. Second half round iteration: from check node and 
initial value to symbol node Cc{n) 



'J^ —(On — l) k 

nm ~ Pk ^' 



m'eN^\^ 



tnm — {^nm^ ^ ^ ^q') — ( (1/ ^nm)^nm^ ^ ^ ^q' ) 5 

(20) 

where p/e denotes the average probability of occurrence of 
element k and On indicates column weight for the n-th symbol 
node. Let Pn = {pk^ k G Zq^}, e.g., we have = {|, ^, |} 
for A^n^ mapping. Note that in ([2Qb the firstly mentioned ix^^ 
denotes the corresponding value generated in Step 1. Here 
it operates as extra information and takes part in the Tanner 
graph. 

4. Soft decision 

Un=Pr^^nm U ^mn] 
meNn 

(21) 

Here also denotes the corresponding value generated in 
Step 1. 

5. Hard decision 



= arg max 

M':{k'} 



(22) 



M:{k}M'-{k'} 



If Cc satisfies the applied mapping's check-relation- 
tab for virtual encoder or the number of iterations exceeds a 
certain value, then the algorithm stops, otherwise we go to 
Step 2. Note that M' mapping denotes the second mapping 
used in hard decision, which may be different from the Ai 
mapping applied in soft iterations. Fortunately, the cardinality 
of the mapping is always far less than that of the M. 

^Note that the definition of A4ni niapping, which will be used as an exam- 
ple in this subsection, is similar to the definition described in Subsection lV-Al I 



mapping. That is to say, the satisfaction of the check-relation- 
tabs will be checked more effectively. For Mni mapping, 
M'ni mapping is selected as {b} {0} and {a,c} {1}. 
Namely, {k} = {a, 6, c} and {k^} = {0,1}. Therefore, 
P(0) = P{b) and P(l) = P(a) + P(c) should be carried 
out in hard decision. Therein, P{i) denotes the probability 
of occurrence of the element i, and {0,1} indicate the PLNC 
mapped symbols which will be transmitted in BC phase. 

So far, the whole PCD algorithm for arbitrary PLNC map- 
ping is presented. 

C. Convergence behavior 

There are many methods that can be used to investigate the 
convergence behavior of iterative decoding. Examples are the 
density evolution algorithm ||38 1 and the extrinsic information 
transfer (EXIT) chart [|39|, both of which are suitable for 
LDPC codes over Gaussian channels. However, the virtual 
LDPC code, namely, check-relation-tab in the considered 
TWR model is different from conventional LDPC codes. 
Moreover, the proposed check-relation-tab is determined by 
the selected PLNC mapping. That is to say, several check- 
relation-tabs should be selected adaptively at the relay. To 
the best of our knowledge, there is no any literature to solve 
the analogous problem up to now. Similar to |24|, we resort 
to simulations in Section VII for confirming the convergence 
behavior of the proposed PCD algorithm. 

D. Complexity analysis 

Since decoding complexity of the proposed PCD approach 
is basically determined by the generated check-relation-tabs, 
we focus complexity calculation on the number and size of 
the check-relation-tabs. Note that the proposed check-relation- 
tab is generated by certain two correlative rows with the 
same index of non-zero elements. The two correlative rows 
is defined as a certain two-rows, in which each row has the 
same row index of the respective LDPC code, e.g., two 0*^ 
rows of Ta and Tb- Let M x N, dr and denote matrix 
size, maximum row weight and the hi-th row weight in degree 
distributions of an arbitrary irregular LDPC code, respectively. 
The number of the check-relation- tabs for a given PLNC 
mapping, denoted as Nt, is upper bounded by 



Nr 



(23) 



where Xr^ is a ratio of the number of rows with row weight 
to M. The equality can be reached when the non-zero elements 
of both two correlative rows are different from each other. 

We also obtain the size of the check-relation-tab for virtual 
encoder with a given PLNC mapping as 



Se 



(24) 



where q' denotes the range of PLNC mapped symbols at the 
relay, q < q' < q^. 

Furthermore, the size range of the check-relation-tab for 
PCD decoder is obtained as 



q < OD <q 



(25) 
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The number and size of the check-relation-tabs are de- 
termined by the selected PLNC mapping and the applied 
optimizations. We will state it in Section [Vll in detail. 

VI. Two-stage CNC mapping with PCD 

Recall that the design criterion of M{Ca^Cb) is to 
maximize the MED, as stated in Section IV. However, it 
is difficult to directly maximize the MED. In this section, 
we shall propose a two-stage CNC mapping based on PCD 
decoding, which optimizes the symbol distance first and then 
the Hamming distance. 

A. Two-stage CNC mapping with PCD (TS-CNC-PCD) 

In this method, the traditional CNC mapping proposed 
in ||25| is divided to two steps by a maximum splitting 
and minimum merging (MSMM) strategy for maximizing 
the symbol distance first and then a minimum correlation 
optimization (MCO) method based on the proposed PCD 
approach is presented to optimize the Hamming distance. 

1 ) Channel coding structure at the relay: The proposed TS- 
CNC-PCD is composed of the first CNC mapper (GF((7+p2)), 
PCD decoder (GF((7+p2)), the second CNC mapper (GF((7+ 
Pi)) and {q +pi)PSK/QAM modulator as shown in Fig. [3l 
Upon receiving the superimposed signal, the relay initializes 
the soft value for PCD decoder by taking into account the 
1st CNC mapping. The output of the PCD decoder, i.e., the 
2nd CNC mapped codeword packet, Cc = A1(C^,Cb), is 
then modulated to {q +pi)PSK/QAM accordingly to obtain 
Xc. The cardinality of the received symbol pairs at the relay 
node is q^ . With the first application of the PLNC which is 
denoted as the 1st CNC mapping, we reduce the cardinality 
to q -\- P2 by clustering some received symbol pairs into one 
decoding symbol. Next, we make further reduction in the 
cardinality through the second application of the PLNC, i.e., 
the 2nd CNC mapping. Wherein, we cluster some decoding 
symbols into one broadcasted symbol. Here, both pi and 
P2, Pi^P2 ^ ^q'^-q^ represent the possible expanding of the 
cardinality compared to q. For example, we have pi = 0(1), 
P2 = 5(8) for the QPSK(5QAM) at the relay node when 
QPSK (q = 4) is used at two sources. Note that the 1st and 
2nd CNC mappings denote the symbol pair mappings while 
the aforementioned M mapping indicates the codeword pair 
mapping. 

2 ) Two-stage CNC mapping: Since the best CNC mappings 
for any choice of constellation size have been presented 
in ll25ll . we just use for reference and generate the two-stage 
CNC mapping. Replacing x Z4 by Zq x Zq, we can easily 
obtain the best CNC mappings for g-ary modulations using 
the Algorithm 1 in |25|. We enlarge the MED between two 
distinct codewords (e.g. and C[^) by extending the symbol 
distance between two individual coded symbols (namely CNC 
mapped clusters), if the codewords are considered instead of 
uncoded symbols. Using the proposed MSMM strategy, we 
divide the traditional CNC mapping (denoted as M^) into 
the 1st CNC mapping (denoted sls M^) and the 2nd CNC 
mapping (denoted as A^^) for soft value initialization and hard 
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(a) 1st CNC Mapping (Soft Value Initialization) 
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(b) 2nd CNC Mapping (Hard Value Decision) 

Fig. 4: Two- stage CNC mapping according to the channel ratio 

Hbc/Hac = 7(cos^ -\- j sinO) when q = 2^. 



value decision, respectively. For certain Hac and Hbc, the 
proposed MSMM strategy is formulated as 

a) Maximum splitting (A4^ A4^): Any two symbol pairs, 
belonged to slM^ mapped cluster, should be split into two 
mapped clusters, if the distance between which is greater than 
dmax- All mapped clusters, split from an identical 
mapping, generate sl mapping. Note that two distinct 
mappings can generate an identical mapping for the same 
channel condition, while an identical mapping can be 
divided into more than one distinct mappings for different 
channel conditions. 

b) Minimum merging (M^ M^): Any two distinct 
mapped clusters, split from a mapped cluster, should be 
merged into sl mapped cluster. All mapped clusters, 
merged by an identical mapping, generate sl mapping. 
Any two distinct mappings also should be merged into 
one mapping if the minimum distances of which are both 
equal to dmin- 

Therein, dmax and dmin denote the maximum and minimum 
values among whole distances between any two mapped 
clusters, respectively. A mapping is composed by several 
clusters. Several symbol pairs form a cluster. 

For instance, we consider the 4-ary LDPC codes with the 
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Fig. 3: Partial decoding model at the relay node for TWR fading channels. 



TABLE II: Two-stage CNC mappings for 4-ary LDPC codes with QPSK modulation 



(0, 0) (0, 1) (0, 2) (0, 3) 



(1,0) (1,1) (1,2) (1,3) 



(2, 0) (2, 1) (2, 2) (2, 3) 



(3, 0) (3, 1) (3, 2) (3, 3) 



QPSK modulation. The corresponding best CNC mappings 
have been presented in the Table I in |25|. As depicted in 
Fig. |4] and Table HIl we generate 12 1st CNC mappings 
and 6 2nd CNC mappings, denoted as A^|,i G [0,11] and 
Mjjj G [0,5]. Take Fig. [5]-Case I for example, we have 
four clusters according to the Table I in 1251 . e.g., symbol 
pairs {(0, 1), (1, 0), (2, 3), (3, 2)} should be clustered together. 
However, in two- stage CNC mapping, we only group together 
some symbol pairs (e.g. (0,1) and (1,0) in Mq) due to that the 
distance between which is much smaller than dmax- Whereas, 
we classify the other symbol pairs (e.g. (2,3) and (3,2)), which 
are more far away from (0,1) or (1,0), into another independent 
cluster. It is because that the distance between these two 
clusters {(2,3),(3,2)} and {(0,1),(1,0)} is greater than dmax- 
Certainly, these separated clusters (e.g. (0,1), (1,0), (2, 3) and 
(3,2)) are merged to one cluster again in Mq. Similarly, 
MpMi and MyM2 are generated according to the pro- 
posed MSMM strategy as depicted in Fig. [5]-Cases II and III. 
Here, to avoid confusion, we let {a',b\c\ . . .} indicate the 
broadcasted symbols {0,1,2,...}, which are not same as the 
decoding symbols {a, 6, c, . . . } during the period of operating 
PCD algorithm. 



3) Check-relation-tabs and TS-CNC-PCD: Note that only 
the 1st CNC mappings are operated before soft iteration in 
the proposed PCD approach, we just need to generate the 
check-relation-tabs according to A^^. However, the additional 
check-relation-tabs for also need to be generated if 

the satisfaction of the check-relation- tabs should be checked 
during each iteration. Similar to Subsection V-A and V-B, we 
can obtain the check-relation- tabs for all 1st and 2nd CNC 
mappings easily and operate the proposed PCD algorithm 
directly, although we should modify the cardinality of the 
designed PLNC mapping here compared with that of Section 
V. 

For instance, we have Ml,i G [0, 11] if the 4-ary LDPC 
codes and the QPSK modulation are applied at two source 
nodes. Take the M% mapping of Table [III for example, 
which is displayed in Fig. [5]-Case III. Each symbol pair 
(Ca(^), is mapped to one of 12 elements based on 

the fading conditions Hbc/Hac — (l+j)/2. All the possible 
mappings are listed below: 

{(0, 1), (1, 2)} ^ {a}, {(0, 3), (2, 0)} ^ {6}, 

{(1, 3), (3, 0)} ^ {c}, {(2, 1), (3, 2)} ^ {d}, 
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Case I: H.jH,^ = 1 Case II: H,jH^, = j Case III: H.jH^^ = (l + j)/2 

Fig. 5: Received signal constellation with two-stage CNC mapping at the relay. 



{(0, 0)} ^ {e}, {(0, 2)} ^ {/}, {(1, 0)} ^ M, 

{(1, 1)} ^ W, {(2, 2)} ^ {i}, {(2, 3)} ^ {j}, 

{(3,1)}^ W, {(3,3)}^ W. 

Each element constrained by the check functions corre- 
sponds to one symbol pair in ^ ^ M^}, 
the range of which is in {a, 6, c, e, /, z, j, /c, Accord- 
ing to the 2nd CNC mapping ^ we should operate P(a') = 
P{f) + ^(^) + PU) + ^(^), ^(^0 = P{e) + P(d), P(cO = 
P{a) + P(0, = P{b) + P(/i), P(eO = P{c) + ^(^) in 
hard decision. 

B. Minimum correlation optimization (MCO) 

Due to the irregular cardinality of the 1st CNC mappings, 
the size of mostly generated check-relation- tabs for PCD 
decoder can approach the maximum value, i.e., Sd = q'^'^- 
That is to say, we can not obtain larger coding gains from 
these check-relation-tabs. Note that the check-relation-tabs 
are fixed once a PLNC mapping and the correlative rows of 
two LDPC codes (Ta and Tb) are all selected. Therefore, a 
kind of correlative rows optimization, named as the minimum 
correlation optimization (MCO), is proposed for selecting the 
better check-relation-tabs. 

Using the proposed MCO method (Algorithm [TJ, we obtain 
a collection of non-zero elements distributions for certain cor- 
relative rows of and Tb- Among these distributions, some 
special correlative rows like {t^,...,//},/^ G Zq, can generate 
the check-relation-tabs with a relatively small size. A check- 
relation-tab with a smaller size always has a lager Hamming 
distance. Certainly, larger coding gains may be obtained by 
these check-relation-tabs. Moreover, these special correlative 
rows can also reduce the number of the check-relation- tabs 
sharply. In theory, it is likely to generate different check- 
relation-tabs when the row weight or the non-zero elements 
of any two rows of each LDPC code are not similar to each 
other. However, only one check-relation-tab for the virtual 
encoder and the corresponding one check-relation-tab for PCD 
decoder needs to be obtained for one or more selected adaptive 
PLNC mapping when the regular LDPC codes are applied 
and the non-zero elements in every row follows some special 
pattens like {7^,.. .,7^},/^ G Zq. Considering the trade-off of 
complexity and performance, we select the special correlative 
rows like {7^, . . . , 7^}, 7/ G Zq. 



Data: given r^, q, q\ M 

Result: a collection of non-zero elements distribution for 
the correlative rows 

1 Initialization: {r]i,r]2, . . . , 7^^^} = {1, 1, . . . , 1}, 

{6,6,...,^rJ = {i,i,...,i} ; 

2 Set the expected collection: Cexp = 0; 

3 Set the maximum average number of possible generation: 

C-max — Q -> 

4 while 7^1 < g do 

5 G M, generate the relevant check-relation tabs; 

6 Let the temporary value Gmaxtemp ^Qual to the 
maximum average number of possible generations of 
the lastly check-relation tabs; 

7 if Gmaxtemp — Graax then 



8 




if Gmaxtemp — Gmax then 


9 






Add {7?i,7?2,...,^rJ,{6,6,---,^rJ ^ 

Gexp 1 


10 




else 


11 






Set Gfnax — Gqjiaxtemp^ G^xp — 0? add 


12 




end 


13 


end 




14 


if ^n(^n) = q, ne [2,r^] then 


15 




iniVn) = 1, ^n-l(^n-l) = ^n-l(^n-l) + 1; 


16 


else 




17 




^r. = ^r. + 1; 


18 


end 




19 


if 6 


= q then 


20 




{6,6,---,^rJ = {1,1,...,!}, Vr^ = Vr. + 


21 


end 




22 end 







Algorithm 1: Minimum correlation optimization (MCO) 



C. Simplification of TS-CNC-PCD 

To simplify the decoding process, the two-stage CNC map- 
ping can be integrated into one- stage mapping. Specifically, 
the identical mapping is used for soft value initialization 
and hard value decision. As shown in Fig. [3l the simplifi- 
cation of TS-CNC-PCD is composed of an adaptive CNC 
mapper (GF{q + pi)), PCD decoder (GF{q + pi)) and 
{q +Pi)PSK/QAM modulator. According to the same design 
principle, we can generate the corresponding check-relation- 
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tabs. The expected error performance may be worse. However, 
the complexity (size of check-relation-tab) is decreased dra- 
matically. More details can be found in |JJ, L2J. 

VII. Simulation results 

In this section, we present some simulation results to 
illustrate the convergence behaviors and the coding gains of 
the proposed PCD approach. For simplicity, each node uses 
the same transmission power of one unit and observes the 
same noise power given by a^. Define an average SNR per 
information symbol as 2^^, where r is the channel code rate. 
According to the proposed MCO method, we generate a 4- 
ary code from a binary LDPC code "252.252.3.252", which 
is produced by Mackay [40 1, through replacing {1, . . . , 1} by 
{t^,...,/^}, t] e Z4. Code length, code rate, row weight and 
column weight are 504, 0.5, 6 and 3, respectively. 

The proposed partial decoding method, TS-CNC-PCD, is 
simulated in kinds of TWR channels. The selections for the 
two- stage CNC mapping and traditional CNC mapping are 
based on instantaneous realizations of the channel gain pairs 
{Hac^ Hbc} using Fig. |4] (a,b) and Fig. |4] (b), respectively. 

For comparison, two benchmark systems are considered. 
One is the uncoded case, where QPSK modulation is applied 
and the relay demodulates using different types of PLNC 
mappings. The other is the coded conventional XOR case, 
where the same 4-ary LDPC is applied at the source and 
the relay performs traditional BP decoding with conventional 
XOR mapping. 

A. Deterministic channels 

Fig. [6] shows the error performance of the proposed PCD 
algorithm over three deterministic channels whose channel 
coefficients are fixed. The SER plotted in simulation is defined 
at the relay node over the MA phase only. The number of 
maximum iterations is set as 10, 20 and 30. In general, 
larger maximum iteration leads to better performance for the 
proposed TS-CNC-PCD method no matter which channel gain 
pair is considered, as shown in Fig. [6l 

In Fig. Oa) where Hac = Hbc = 1, it is observed 
that uncoded XOR and CNC have the same performance. It 
is because the CNC mapping is degraded to XOR mapping 
when the two channel gains Hac and Hbc are identical. 
Interestingly, the proposed TS-CNC-PCD is a little bit better 
than the XOR-BP. Since the traditional XOR can work well 
in this case, this observation confirms that there is no inherent 
loss coupling this non-linear operation with a LDPC code 
(i.e. TS-CNC-PCD) compared to the referred linear operation 
XOR-BP. Due to that the non-linear operation usually result 
in shattered pairwise check constraints. In other words, it 
decreases the Hamming distance of generated codeword space. 
Fortunately, the proposed TS-CNC-PCD not only compensates 
the aforesaid performance loss but also obtains extra im- 
provement, through enlarging the cardinality of the decoding 
symbols and increasing the symbol distance between two 
distinct PLNC mapped symbols. It is clear that the coding 
gains of all considered coded systems are more than 3.5dB 
compared to the uncoded scenarios at SER = 10~^. 
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Fig. 6: Convergence behaviors of the proposed PCD algorithm. 
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Fig. 7: Performance comparisons in the TWR block fading 
channels. 

From Fig. [6tb), we can see that the conventional XOR 
mapping does not work anymore when Hac = 1 and 
Hbc = j- The reason is that in this case the XOR mapping 
badly decreases the symbol distance between two distinct 
PLNC mapped symbols although it has no loss in Hamming 
distance of generated codeword space. Fortunately, the coding 
gains of TS-CNC-PCD is also more than 3.5dB compared to 
the uncoded CNC at SER = 10"^ 

Lastly, as shown in Fig. Oc) where Hac = 1 and 
Hbc = (1 + j)/2, the XOR-BP coding method outperforms 
the uncoded CNC more than 1.25dB at SER = 10"^ while 
the uncoded XOR mapping still does not work well. At the 
same time, the proposed TS-CNC-PCD obtains more than 2dB 
coding gain compared to the XOR-BP at SER = 10~^. 

In general, our proposed TS-CNC-PCD not only obtains the 
best performance but also achieves the constant coding gain 
in considered TWR deterministic channels. 

B. Block fading channels 

Suppose that the channel gains on all links follow Rayleigh 
distribution and are independent. We assume £^[|i^ACp] = 
E[\Hbc\'^] = 1, where notation E[-] denotes expectation func- 
tion. The black solid lines, denoted as "Outage Probability", 
are actually the lower bound of the outage probability of the 
TWR block fading channels according to ([5]). The maximum 
number of iterations is fixed at 30 for the coded cases. 

Fig. [7] shows the frame error rate (FER) performance of 
the Rayleigh channels. For the uncoded cases, the CNC 
outperforms the XOR about 4 dB at FER = 2 x 10"^. At 
the same FER, the coding gain of the coded XOR is about 
8.5 dB. Moreover, the coding gain of the coded CNC is about 
6 dB at FER = 2 x 10~^. The reason that the coding gain of 
TS-CNC-PCD is less than that of XOR-BP is that the uncoded 
denoising mapping used at the relay node can eliminate part of 
the noise, as confirmed in ll25ll . Nevertheless the TS-CNC-PCD 
still outperforms the XOR-BP about 2 dB at FER =2 x 10"^. 
At the same time, the gap between the TS-CNC-PCD and 
the lower bound of the outage probability is 4.5 dB. By any 



possibility, we can reduce the gaps by increasing the code 
length. 

Lastly, it is important to see that both the relay and the 
sources obtain the same FER performance. These phenomena 
confirm that the errors at the relay have severe impact on the 
performance of whole system. This further confirms the needs 
for advantage decoding at the relay, such as the proposed PCD 
algorithm. 

VIIL Conclusion 

In this paper, maintaining the traditional channel coding 
structure at the two sources, we propose a general relay 
decoding framework, called pairwise check decoding, or PCD, 
for any given PLNC mapping. The check-relation-tab for the 
superimposed LDPC-coded packet pair at relay is formed 
according to the selected PLNC mapping. Our proposed PCD 
algorithm is universal for any adaptive PLNC mapping. In 
order to optimize the MED, we also present a partial de- 
coding method at the relay, TS-CNC-PCD, for LDPC coded 
TWR block fading channels. Moreover, a check-relation-tab 
optimization MCO is introduced to improve performance. 
Simulation results confirm that the proposed TS-CNC-PCD 
has significant coding gains compared to the conventional 
XOR with BP and the uncoded system for certain TWR 
deterministic channels. For TWR fading channels, the TS- 
CNC-PCD also considerably outperforms the conventional 
XOR with BR 
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